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Aims. The properties of the broad-band emission from the high-frequency peaked BLLac H 2356— 309 (z = 0.165) are investigated. 
Methods. Very High Energy (VHE; E > 100 GeV) observations of H 2356—309 were performed with the High Energy Stereoscopic 
System (HESS) from 2004 through 2007. Simultaneous optical/UV and X-ray observations were made with the XMM-Newton 
satellite on June 12/13 and June 14/15, 2005. NRT radio observations were also contemporaneously performed in 2005. ATOM 
optical monitoring observations were also made in 2007. 

Results. A strong VHE signal, ~13<r total, was detected by HESS after the four years HESS observations (116.8 hrs live time). The 
integral flux above 240 GeV is I(>240 GeV) = (3.06 ± 0.26 st at ± 0.61 syst ) x 10~ 12 cm^ 2 s"\ corresponding to ~1.6% of the flux 
observed from the Crab Nebula. A time-averaged energy spectrum is measured from 200 GeV to 2 TeV and is characterized by a 
power law (photon index of T= 3.06 ± 0.15 s t a t ±0.10 sys t). Significant small-amplitude variations in the VHE flux from H 2356— 309 
are seen on time scales of months and years, but not on shorter time scales. No evidence for any variations in the VHE spectral 
slope are found within these data. The XMM-Newton X-ray measurements show a historically low X-ray state, characterized by 
a hard, broken-power-law spectrum on both nights. 

Conclusions. The broad-band spectral energy distribution (SED) of the blazar can be adequately fit using a simple one-zone 
synchrotron self-Compton (SSC) model. In the SSC scenario, higher VHE fluxes could be expected in the future since the observed 
X-ray flux is at a historically low level. 

Key words. Galaxies: active - BL Lacertae objects: Individual: H 2356— 309 - Gamma rays: galaxies 
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1. Introduction 

H 2356-309 



was 



Costamante fc Ghisellini (2002) 
for 



suggested 
as a strong 
emission, based on 



by 

candidate 
simple 



the detection of VHE 
one-zone synchrotron self-Compton (SSC) modeling of its 
spectral energy distribution (SED). Motivated by these 
predictions, the HESS collaboration observed this blazar 
in 2004, resulting in the discovery of VHE 7-ray emission 
(|Aharonian et al. 2006a|l . Although several blazars are 
known to have variable VHE 7-ray fluxes (see, e.g., 
lAharonian et al. 2007aP and photon spectra (see, e.g., 
Aharonia n et al. 2 005 ). no strong evidence was found for 
such variability in the HESS data on this source from June 
through December 2004. 

Hosted by an elliptical galaxy at a redshift of z = 0.165 
(jFalomo 1 991). H 2356— 309 is among the more distant 
blazars detected at VHE energies. As a result, its observed 
spectrum is expected to be strongly affected (softened) by 
the absorption of VHE photons on the extragalactic back- 
ground light (EBL). The unexpectedly hard (r - 3.1) VHE 
spectrum measured by HESS in 2004 was used to confirm 
the strong constraints on the density of the EBL, near 
lower limits based on the integrated light of resolved galax- 
ies ([Aharonia n et al. 2006a[) . given by the HESS spectrum 
from another VHE blazar (1ES 1101-232). 

At X-ray energies , H 2356—309 was initi ally detected by 
the UHURU satellite (Forman et al. 1978)] and later by the 
HEAO-I satellite |(Wood et al. 1984)| It is among the X- 
ray brightest of the 150 high-frequen cy peaked BL Lac ob- 
jects (HBL) in the Sedentary Survey (Giommi et al. 2005)| 
which contains HBLs characterized by their exception- 
ally high X-ray-to-radio flux ratio. BeppoSAX observa- 



(Costamante et al. 2001) of H 2356-309 in June 1998 



tions 

yielded an X-ray spectrum extending up to 50 keV, char 
acterized by a broken power-law shape (ri = OJSIqq^; 
T 2 = 1.10 ± 0.04) with a break at 1.8 ± 0.4 keV. A softer 
power-law spectrum between 2 and 10 keV with an X- 
ray photon index of (Tx = 2.43 ± 0.11) was measured 
(|Aharonian et al. 2006b|) during observations by the RXTE 
satellite on November 11, 2004. The RXTE spectrum is 
softer above 2 keV and the flux (9.7t?'| x 10~ 12 erg cm~ 2 
s _1 ) is ^3 times lower than that measured by BeppoSAX 
(2.5 x 10 -11 erg cm -2 s _1 ) in the same energy band. 

The 2004 RXTE data were taken simultaneously with 
HESS VHE, ROTSE-IIIc optical, and Nancay radio- 
telescope observatio ns. The double-peaked SED derived 
from this campaign (|Aharonian et al. 2006b[) was the first 
to include a contemporaneous VHE detection, and hence 
the first to sample the higher-energy SED peak for this 
object since H 2356-309 was not detected by the EGRET 
satellite. Modeling of these multi-wavelength observations 
demonstrated that the SED of H 2356—309 could be ade- 
quately fit by a simple one-zone stationary SSC scenario. 

Results from 76.9 hrs of new HESS observations of 
H 2356-309 in 2005, 2006 and 2007 are reported here. In 
addition, a re-analysis of the previously published 2004 
HESS data (39.9 hrs), with an improved calibration of 
the absolute energy scale of the detector, is presented. 
This improvement in the HESS calibration was not in- 
cluded in the previously mentioned modeling of the 2004 
multi-wavelength observations, therefore this modeling is 
repeated. In addition, the XMM-Newton satellite simul- 
taneously observed H 2356— 309 on two nights (~5 hours 



each night) at X-ray energies, as well as at optical through 
ultraviolet wavelengths, during the HESS observations in 
June 2005. The results from this multi- wavelength observa- 
tion campaign, which has much improved X-ray and optical 
measurements, are also discussed. 



2. HESS Observations and Analysis Technique 

From 2004 through 2007, a total of 175.3 hours of HESS ob- 
servations were taken on H 2356—309 in wobble mode, where 
the telescope tracking position is offset by ±0.5° in either 
right ascension or declination relative to the source loca- 
tion. After applying the standard HESS data-quality selec- 
tion, the exposure is 116.8 hrs live time at a mean zenith 
angle of Z mcan = 19°. The HESS data are calibrated fol- 



lowing the procedures described in |Aharonian et al. (2004) 
The analysis is performed using the standard HESS tools 
and 7-ray-like events are selected using the standard cuts 
(jBenbow 20051 lAharonian et al. 2006cP . The average post- 
analysis energy thresholcfl is 240 GeV at Z mcan . The sig- 
nal (on-source data) is taken from events with recon- 
structed direction falling within a circular region of ra- 
dius 9 cu t=0-ll° centered on H 2356— 309. The background 
(off-sour ce data) is estimat ed using the Reflected- Region 
method |(Berge et al. 2007)] hence the on-source and off- 
source data are measured simultaneou sly from events in 



the same field of view. Equation 17 in Li & Ma (1983) 



is 



used to calculate the significance of any excess. 

The optical throughput of HESS decreases over time 
due to degradation of the reflective surfaces of the HESS 
mirrors and Winston cones, as well as to the accumula- 
tion of dust on the optical elements. Compared to a newly- 
commissioned HESS telescope (mirror installation between 
October 2001 and August 2003), this decrease is on aver- 
age ~26% for the entire data sample, and is 20%, 27%, 
32% and 35% for the H 2356-309 data in 2004, 2005, 2006 
and 2007, respectively. This decrease causes images from 
gamma rays of equivalent energy to appear less bright in 
later data, requiring a renormalization of the event en- 
ergy. This normalization is made using the ratio of effi- 
ciencies determined from simulated and observed muons 
(| Aharonian et al. 2006cp . and eliminates long-term varia- 
tions in the absolute energy scale of the HESS analysis due 
to a changing optical throughput. 

3. HESS Results 

A total of 1185 excess events (12.6 a) is measured from the 
direction of H 2356— 309 in the complete data set. Figure [1] 
shows the on-source and normalized off-source distributions 
of the sq uare of the an gular difference between the source 
position (Falomo 1991) and the reconstructed shower posi- 
tion (6 2 ) for all observations. The observed signal can be 
seen at small values of 8 2 and the distribution of the excess 
is similar to that expected from a simulated 7-ray point 
source. 

H 2356— 309 is clearly detected with a statistical signif- 
icance of more than 5 standard deviations (cr), in each of 
the first three years (2004-06) of HESS observations, and a 
marginal excess (3.4 a) is observed during the brief observa- 
tions in 2007. Detailed results of the HESS observations of 



1 The threshold is corrected to account for the decreased op- 
tical efficiency of the HESS mirrors over the observation period. 
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Table 1. Results from four years of HESS observations of H 2356— 309. 
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06/2004 


53172 


53186 


7.9 


835 


7778 


0.0938 


106 


3.7 


3.86 ± 0.88 


2.1 


8.9 , 11 


0.63 


07/2004 


53201 


53209 


3.0 


280 


2655 


0.0938 


31 


1.8 


1.19 ± 1.18 


0.6 


14 4 


0.84 


09/2004 


53256 


53269 


6.9 


759 


6015 


0.0937 


195 


7.4 


6.91 ± 0.94 


3.7 


4.0 ,6 


0.67 


10/2004 


53286 


53296 


9.9 


950 


8496 


0.0914 


174 


5.7 


4.33 ± 0.80 


2.3 


9.9 ,6 


0.13 


11/2004 


53316 


53319 


4.1 


300 


2700 


0.0912 


54 


3.2 


3.82 ± 1.22 


2.0 


7.0 , 3 


0.07 


12/2004 


53341 


53353 


8.2 


559 


5867 


0.0926 


16 


0.6 


1.23 ± 1.46 


0.7 


7.2,5 


0.21 


06/2005 


53530 


53541 


20.4 


1200 


12275 


0.0921 


69 


2.0 


1.78 ± 1.04 


1.0 


6.3,9 


0.71 


07/2005 


53556 


53570 


16.5 


1455 


13375 


0.0922 


221 


5.9 


3.08 ± 0.60 


1.6 


12.9 , 13 


0.46 


08/2005 


53583 


53587 


7.3 


495 


4920 


0.0927 


39 


1.7 


0.56 ± 0.85 


0.3 


5.7,4 


0.22 


09/2005 


53614 


53615 


2.6 


152 


1477 


0.0946 


12 


1.0 


2.70 ± 1.67 


1.4 


0.0 , 1 


0.96 


07/2006 


53937 


53942 


8.1 


598 


5926 


0.0893 


69 


2.8 


2.27 ±0.89 


1.2 


2.5,5 


0.77 


08/2006 


53970 


53978 


2.2 


188 


1427 


0.0894 


60 


4.8 


7.39 ± 1.84 


3.9 


0.5,3 


0.91 


09/2006 


53999 


54005 


12.9 


626 


6044 


0.0931 


64 


2.7 


1.26 ±0.77 


0.7 


7.2,6 


0.30 


07/2007 


54295 


54307 


7.0 


502 


4670 


0.0915 


75 


3.4 


3.48 ± 1.07 


1.9 


8.9,10 


0.54 


2004 


53172 


53353 


39.9 


3683 


33511 


0.0927 


576 


9.6 


4.06 ± 0.42 


2.2 


19.1 ,5 


0.002 


2005 


53530 


53615 


46.7 


3302 


32047 


0.0924 


341 


5.9 


2.22 ±0.44 


1.2 


6.1,3 


0.11 


2006 


53937 


54005 


23.2 


1412 


13397 


0.0910 


193 


5.1 


2.36 ±0.56 


1.3 


9.5,2 


0.01 


2007 


54295 


54307 


7.0 


502 


4670 


0.0915 


75 


3.4 


3.48 ±1.07 


1.9 






Total 


53172 


54307 


116.8 


8899 


83625 


0.0923 


1185 


12.6 


3.06 ± 0.26 


1.6 


11.0,3 


0.01 



a The 20% systematic error on the observed integral flux above 240 GeV is not shown. 

The Crab Nebula flux percentage is calculated relative to the HESS value above 240 GeV (|Aharonian et al. 2006c)) . 
c The x 2 , degrees of freedom (NDF), and corresponding \ 2 probability P(x 2 ) are given for fits of a constant to I(>240 GeV) 
binned nightly within a dark period, or monthly within a year, or yearly within the total. 
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Fig. 1. Distribution of 8 2 for on-source events (points) and 
normalized off-source events (shaded) from observations of 
H 2356—309. The dashed curve represents the expected dis- 
tribution for a VHE 7-ray point source with photon index 
T = 3.06 at a zenith angle of 20°. The vertical line de- 
notes the on-source integration region relative to the nom- 
inal source position fF alomo 199 If) . 
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Fig. 2. Time-averaged (2004-2007) VHE photon spectrum 
of H 2356—309. The dashed line represents the best x 2 fit 
of a power law to the data up to 2 TeV, and a linear ex- 
trapolation at higher energies. The upper limits are at the 
99% confidence level (Feldman & Cousins 1998) Only the 
statistical errors are shown. 



H 2356— 309 for individual observation periods are summa- 
rized in TableQ] The table contains the dead-time-corrected 
observation time, the number of on and off-source events, 
the on/off normalization (a), the observed excess and its 
statistical significance for various intervals (e.g. dark peri- 
ods) in which H 2356—309 was observed. 

A two-dimensional Gaussian fit to the sky map of the 
observed excess demonstrates that the excess is point-like 
with an upper limit (99% confidence level) on the ex- 
tension of 1.3'. The centroid of the fit, HESS J2359-306, 
is located at a J20 oo = 23 h 59 m 7.8 s ± 1.6* tat ± 1.3J yst and 
£j2ooo = -30°37'18" ± 20£ at ± 20" t . HESS J2359-306 is 



spatially consistent with the position (jFalomo 1991[) of the 
blazar (a J20 oo = 23 h 59 m 7.8 s , <5 J20 oo = -30°37'38") as ex- 
pected. The distance between the two positions is 20"±31". 

3.1. VHE Spectrum of H 2356- 309 

Figure [2] shows the four-year (2004-2007) time-averaged 
photon spectrum of H 2356—309. The best x 2 fit of a power 
law (dN/dE = I a (E/l TcV)~ r ) to these data yields a pho- 
ton index T = 3.06 ± 0.15 sta t ± 0.10 sys t, and a \ 2 °f 5.7 
for 6 degrees of freedom. The last spectral point at ~1.7 
TeV is not significant (0.8 a), however removing it from 
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Table 2. Results of power-law fits to the VHE spectra of H 2356-309. 
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no - 


13 


cm s 


-l 
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2004 (AH06)^ 


0.165 


1.041 


3.06 


± 


0.21 s tat 


± 


0.10 Syst 


3.08 


± 


0.75 s tat 


± 


0.62 S y S t 


3.9 


6 


0.69 


2004_ 


0.200 


1.500 


2.97 


± 


0.19 st at 


± 


0.10 syst 


4.69 


± 


0.86 stat 


± 


0.94 syst 


7.8 


5 


0.17 


2005~ 


0.200 


2.000 


2.99 


± 


0.39 stat 


± 


0.10 syst 


2.92 


± 


0.89 atat 


± 


0.58 S yst 


5.0 
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0.55 


2006^ 


0.200 


2.000 


3.43 
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0.41 stat 


± 


0.10 syst 


2.13 


± 


0.79 stat 


± 


0.43 syst 


0.7 


2 


0.70 


Total 


0.200 


2.000 


3.06 


± 


0.15 s tat 


± 


0.10 Syst 


3.29 


± 


0.45 stat 


± 


0.66 S y st 


5.7 


6 


0.47 



a The differential flux normalization is calculated at 1 TeV. 

b The AH06 entry is the previously published HESS result (Aharonian ct al. 2006a) for the 2004 data and is not corrected 
for long-term changes in the optical efficiency of the system. 

c This 2004 entry is derived from the exact same data as presented in AH06, but with a correction (see text) applied to 
account for the optical efficiency changes within the data. The 2005, 2006 and total entries also have this correction applied. 

Due to low statistics, the 2006 spectrum is determined using 4 bins/decade in energy instead of the 8 bins/decade used for 
the other photon spectra. 
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Fig. 3. Annual VHE photon spectra observed by HESS 
from H 2356— 309. Each line represents the best \ 2 n t °f 
a power law to the observed data. Only the statistical er- 
rors are shown. 
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Fig. 4. Integral flux, I(>240 GeV), measured by HESS from 
H 2356— 309 during each dark period of observations. For 
each point the time- averaged T = 3.06 is assumed and only 
the statistical errors are shown. The horizontal line repre- 
sents the average flux for all the HESS observations. The 
three horizontal line segments are the average 2004, 2005 
and 2006 fluxes. 



the fit does not appreciably change the result. Fitting the 
observed photon spectrum with a more complex function, 
such as a power-law with a cutoff or break, does not yield 
a \ 2 that is significantly improved (as determined by an 
F-test). 

The high photon statistics observed in 2004, 2005 and 
2006 enable the determination of a spectrum for each of 
those years. Unfortunately a spectrum could not be gen- 
erated for the 2007 data due to low statistics. The annual 
spectra for the first three years of observations are shown 
in Figure [3] together with the best \ 2 n t °f a power law to 
the data. Table [2] gives the fit results for each of the annual 
spectra, as well as for the total spectrum. The epoch, lower 
and upper energy bounds, photon index (r), differential 
flux normalization at 1 TeV (I ), X 2 , degrees of freedom 
(NDF), and \ 2 probability (P(x 2 )) for each fit are shown. 
Although the differential flux normalization I Q is ~1.5 times 
larger in 2004 than in 2005 and 2006, there are no significant 
changes in the spectral slope during the HESS observations. 



3.2. VHE Flux from H 2356-309 

The observed integral fiusQ above 240 GeV for the 
entire data set is I(>240 GeV) = (3.06 ± 0.26 sta t ± 
0.61 syst ) x 10~ 12 cm _2 s _1 , corresponding to ~1.6% of 
I(>240 GeV) determine d by HESS from the Crab Nebula 
(| Aharonian et al. 2006c[) . Figures 0] and [5] show the flux 
measured for each dark period and night, respectively. The 
integral flux I(>240 GeV) for each year of observations, as 
well as for each dark period, is shown in Table Q] Also given 
are the x 2 and corresponding probability for a fit of a con- 
stant to the data when binned by nights within each dark 
period, by dark periods within a year, and by year within 
the total observations. There are clear indications that the 
VHE flux from H 2356—309 varies weakly on time scales of 
months and years. There is no evidence for flux variability 
on any shorter time scale within the HESS data, however 
the average flux is too low to significantly detect similar 
weak variations during such brief exposures. 

2 All HESS fluxes (i.e. annual, dark period and nightly values) 
in this article are calculated using the time-averaged T — 3.06. 
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Fig. 5. Integral flux, I(>240 GeV), measured by HESS from H 2356—309 during each night of observations in 2004 (top), 
2005 (top-middle), 2006 (bottom-middle) and 2007 (bottom). For each point the time-averaged T = 3.06 is assumed, and 
only the statistical errors are shown. The horizontal line in each figure represents the average flux measured by HESS in 
the respective year. The night of the RXTE observations (MJD 53320) is denoted by the dotted vertical line in the 2004 
light curve. The nights (MJD 53534 & 53536) of the XMM-Newton observations are marked by dotted vertical lines in 
the 2005 figure. 
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Fig. 6. Renormalized VHE spectrum measured by HESS 
from H 2356—309 in 2004 compared to the previously pub- 
lished version (lAharonian e t al. 2006a). Only the statistical 
errors are shown. 

3.3. Optical Efficiency Correction 

The previously published (|Aharonian et al. 2006a|) spec- 
trum of H 2356— 309 in 2004 was not corrected for opti- 
cal sensitivity changes. Figure [5] illustrates the effect of 
correcting the energy of individual events for the relative 
optical efficiency of the system. The results of the best fit to 
each set of points is shown in Table [5J As can be seen, the 
corrected spectrum has a significantly larger flux normal- 
ization (/ ), but the photon index (r) is unchanged. The 
corrected 2004 integral flux is ^50% higher than previously 
published. The increase in integral flux due to the optical 
efficiency correction is relatively large for H 2356—309 be- 
cause of the softness of the observed VHE specrum. The 
correction is smaller for harder spectrum sources (e.g., typ- 
ical Galactic VHE sources). 

4. Multi-wavelength Observations 

The emission from blazars is known to span over twenty or- 
ders of magnitude in energy and is variable at all energies. 
Accurate modeling of the underlying processes in these ob- 
jects thus requires simultaneous observations at many wave- 
lengths. Two such simultaneous multi-wavelength observa- 
tion campaigns were organized by the HESS collaboration 
for H 2356-309. The first of these included the RXTE X-ray 
satellite, the ROTSE-IIIc optical telescope and the Nancay 
Radio Telescope (NRT), and is henceforth referred to as 
the RXTE campaign. The second campaign involved X- 
ray and optical/UV observations with XMM- Newton satel- 
lite and further NRT data, and is henceforth referred to as 
the XMM-Newton campaign. In addition, results of ATOM 
optical monitoring of H 2356—309 beginning in November 
2006 are presented. 

4.1. RXTE Campaign 

Results from the RXTE campaign were first reported in 
Aharonian et al. (2006b). As only the normalization of 
the HESS observations have changed, the previously re- 
ported lower-frequency results are summarized here. The 
RXTE/PCA observed H 2356-309 for a total of 2.14 ks on 
November 11, 2004 (MJD 53320) as a target-of-opportunity 
request. The constant X-ray flux was 9.7±^3 x 10~ 12 erg 
cm~ 2 s _1 in the 2—10 keV band. The X-ray spectrum was 



generated by combining data from three separate PCUs, 
yielding a single-PCU equivalent exposure of 5.42 ks, and 
is compatible with a power-law function with photon index 
T x = 2.43 ± 0.11. The average NRT 11-cm flux between 
June 11 and October 10, 2004 was 40 ± 8 mJy. During the 
observations by HESS in 2004, ROTSE-IIIc measured ap- 
parent R-band magnitudes m(R) between 16.1 and 16.9, 
with an estimated contribution from the resolved elliptical 
host galaxy of m{R) = 17. Due to poor weather, HESS was 
unable to observe the blazar on MJD 53320. However, there 
is no evidence for variability in the November 2004 VHE 
flux and this flux is consistent with the average 2004 VHE 
flux. As there is no evidence for changes in the spectral 
slope of H 2356—309 on any time scale, the time-averaged 
VHE spectrum from 2004 is assumed to represent the VHE 
state of the blazar during the RXTE campaign for the SED 
modeling in Section [5J The same assumption was made in 
Aharonian et al. (2006b). 

4.2. XMM-Newton Campaign 

The XMM-Newton satellite was pointed at H 2356-309 on 
two nights in 2005. The first XMM-Newton pointing began 
at 23:59 UTC on June 12, 2005 and ended at 5:17 UTC on 
June 13, 2005 (MJD 53334). The second pointing started 
at 23:57 UTC on June 14, 2005 and finished at 5:19 UTC 
on June 15, 2005 (MJD 53536). Simultaneous HESS obser- 
vations of the blazar were also made during both of the 
XMM-Newton pointings. At radio wavelengths, the 11-cm 
flux from H 2356— 309 was contemporaneously monitored 
by the NRT. 

4.2.1. HESS Results 

The total HESS exposure during the XMM-Newton cam- 
paign is 4.68 hrs live time after selection criteria. The data 
divided approximately evenly over the two nights with ex- 
posures of 2.44 hrs and 2.24 hrs on MJD 53534 and 53536, 
respectively. HESS did not detect a significant excess from 
H 2356-309 on either night of the XMM-Newton observa- 
tions. During these two nights, a total of 302 on-source 
events and 2907 off-source events were measured with an 
on-off normalization of 0.0930, corresponding to an excess 
of 32 events (1.8cr). Due to the lack of a significant detection 
it is not possible to produce a photon spectrum from the 
data. Assuming the observed excess is from 7-rays emitted 
by H 2356-309, the average flux during the XMM-Newton 
observations, I(>240 GeV) = (5.7±2.9 sta t±l.l syst ) x 10~ 12 
cm~ 2 s , is consistent with time-averaged value from 2005. 
The flux for each of the nights (MJD 53534 and 53536) in 
the XMM-Newton epoch can be seen in Figure [5] There 
are no significant variations of the nightly flux, or of the 
run-wise (~28 min) flux in either of the individual nights, 
during the XMM-Newton observations. Since there are no 
indications of any VHE flux variations in 2005, the time- 
averaged VHE spectrum from 2005 is used for the SED 
modeling of the each night of the XMM-Newton campaign 
in Section [5] 

4.2.2. XMM-Newton X-ray Results 

During the two pointings, the XMM-Newton EPIC instru- 
ments were set in timing (PN and MOS1 cameras) and large 
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Instrument 


MJD 


Model 




-f^brcak 

[keV] 


r 2 


Fq. 1-2.0 kcV 

[10~ 12 ergs 


-F2-10 kcV 

cm" 2 s~ x t 


Xr , NDF^ 


BeppoSAX 
RXTE 


51350-51 
53320 


BPL 
PL 


0.78t8;g? 
2.43 ± 0.25 


1.8 ±0.4 


1.10 ±0.04 




25 
9.70i?j 




XMM-Newton^ 


53534 


PL 
BPL 


2.23 ±0.01 
2.08 ±0.03 


1.00 ± 0.08 


2.32 ±0.02 


25.9 
23.2 


8.05 
7.46 


1.41 ,541 
0.95,539 


XMM-Newton c 


53536 


PL 
BPL 


2.16 ±0.01 
1.89 ±0.10 


0.66 ±0.14 


2.20 ±0.02 


25.9 
22.3 


9.54 
9.15 


1.52,541 
1.21,539 



a The reduced Xr an< i degrees of freedom (NDF) correspond to fits of either a power-law (PL) function or broken power-law 
(BPL) function to the XMM-Newton X-ray data. 
b The integrated fluxes have the effects of Galactic absorption removed. 
c All errors from the XMM-Newton fits are quoted at the 90% confidence level. 



window (MOS2 camera) modes. On MJD 53534 the total 
live time of the XMM-Newton exposures are Tpn = 16.5 
ksec for the EPIC PN instrument and Tmos2 = 17.6 ksec 
for the EPIC MOS2 detector. The live times are similar, 
Tpn = 16.6 ksec and T MO S2 = 17.7 ksec, on MJD 53536. 
The XMM-Newton Science Analysis System (SAS version 
7.1.0) is used to process the data with the July 2008 cal- 
ibration files. The spectral and timing analysis are per- 
formed with XSPEC vll.3.2ag and FT00LS V6.3.2. The 
MOS2 data are taken in Large Window mode, while the 
MOS1 in Timing mode. Data from the MOS1 detector are 
not included, as they are in qualitative agreement with the 
other instruments but exhibit much stronger noise levels. 
The 0.1-1 keV and 4-10 keV count rates of the PN, as 
well as the corresponding hardness ratio, are all constant 
in time on each of the two nights. Therefore only the night- 
averaged energy spectra are presented. 

The effects of pile-up in the EPIC instruments were ex- 
plored using epatplot. No significant pile- up is present. 
The spectrum is obtained by selecting only single pixel 
(PATTERN=0) events for the MOS2 data, and single plus 
double pixel events (PATTERN<4) for the PN data. 

For the MOS2 analysis, the signal is extracted from a 
circle of radius 45" centered on the source centroid. The 
background is taken from an annulus around the source 
with an inner radius of 120" and an outer radius set by 
the border of the CCD window. For the PN analysis, signal 
photons are selected from rows 28 < RAWX < 48 and the 
background is estimated from 2 < RAWX < 18. The energy 
range of the PN is restricted to 0.5 to 10 keV, while MOS2 
events between 0.15 and 10 keV are accepted. For the spec- 
tral determination, different re-binning schemes are used, 
always requiring at least 50 counts per new bin. Ancillary 
and response files are produced with rmf gen and arf gen. 

The spectra of the PN and MOS2 detectors are com- 
bined in the following fitting, with a free constant (that re- 
mained within a few percent of unity) to allow for the differ- 
ent normalization between the two instruments. The spec- 
tra are fit with source models including Galactic absorption 
along the l ine of sight to H 2356 —309. The absorption model 
TBab s by Wilms et al. (2000) is used, with cross-sections 
by Verner et al. (1996)| Theabsorbing column density is 
kept fixed at the Galactic value of Ah = 1.44 x 10 20 cm -2 
(jKalberla et al. 2005[) . 



A single power-law model does not provide an accept- 
able fit to the X-ray data of either night (see Table [3]). 
Different models of the Galactic absorption (wabs, phabs) 
were also tested, but none of them provide an accept- 
able fit for a single power-law source model. A broken 
power-law provides an improved fit on both nights (F-test 
> 99.999%), and also with respect to a fit with free A^h 
(reduced xl = 0-95 for 539 NDF vs. xl = L12 for 540 
NDF). 

The results of the broken power-law fits for both nights 
are shown in Table [31 along with a summary of archival X- 
ray spectrum measurements. The measured X-ray spectra 
for the XMM-Newton pointings on MJD 53534 and 53536 
are shown in Figure [7] 

4.2.3. XMM-Newton Optical Monitor Results 

During the two pointings the XMM-Newton optical moni- 
tor (OM) took four exposures with four different filters (V, 
B, U, UVW1). For each exposure, the window on the tar- 
get was set in fast mode. The XMM-Newton OM data are 
processed with xmmsas 7.1.0 using data from the imaging 
mode for the photometry, and from the fast timing window 
for the light curves. No variability is found in any of the 
OM exposures for each of the different filters. The point- 
source analysis of the OM photometry program omsource 
is used to extract count rates, which are converted into 
fluxes following the standard procedures. For the V, B, and 
U filters, the on-source data are taken from an aperture of 
R = 6" and the background is estimated from an annulus 
with 10" < R < 15". An on-source aperture of R = 17". 5 
and background annulus of 20" < R < 25" are used for the 
UVW1 filter. Consistent results were also obtained using 
different background regions. 

The optical spectrum for each of the two nights is 
generated by combining the results from all four filters. 
Each of the V and B band exposures of H 2356— 309 con- 
tains a significant contribution from the host galaxy that 
must be subtracted. The host galaxy is resolved at op- 
tical (Falomo 199U Scarpa et al. 2000 ) and near-infrared 
( |Cheung et al. 2003| ) wavelengths, and is a normal ellip- 
tical galaxy with an effective radius of about 1".7 in the 
V and R bands, somewhat lower in the IR. Using a stan- 
dard de Vaucouleurs radial profile, more than 60% of the 
host-galaxy flux is thus estimated to be contained in the V 
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Fig. 7. X-ray spectra measured during the XMM- Newton observations in 2005. The lefthand and righthand figures 
correspond to the measurements on MJD 53534 and 53536, respectively. The red points (0.5 to 10 keV) and black points 
(0.15 to 10 keV) are data from the PN and MOS2 detectors, repectively. The corresponding fits to the data are also 
shown as lines, and the fit residuals are shown in the lower panels. 



Table 4. XMM -Newton Optical Monitor results. 



Band 




Fi 


T 2 a 


F 2 




[ksec] 


[mJy] 


[ksec] 


[mJy] 


V 


4.3 


0.33 


3.8 


0.32 


B 


3.9 


0.22 


4.3 


0.21 


U 


3.9 


0.24 


4.3 


0.23 


UVW1 


4.4 


0.24 


4.4 


0.23 



a The OM exposures (T) and fluxes (F) are reported for each 
filter on MJD 53534 (Ti & Fi) and 53536 (T 2 & F 2 ). 



4.3. ATOM Results 

Optical observations were taken using the ATOM telescope 



(Hauser et al. 2004)| at the HESS site from November 2006 
on. Absolute flux values are calculated using differential 
phot ometry against several s tars listed in the NOMAD cat- 
alog (Zacharias et al. 2005) A 4 arcsec radius aperture is 
used lor all filter bands. A total of ^500 measurements in 
4 filter bands were taken between November 2006 and the 
end of 2007. The measured apparent magnitudes vary in 
the R band around (16.4 ± 0.2) mag. For the other filter 
bands, the values are (17.1 ± 0.2) mag in B, (17.0 ± 0.2) 
mag in V and (16.1 ± 0.2) mag in I band. 



and B band signal apertures. The host-galaxy contribution 
to the background apertures is negligible. These fractions 
of the host galaxy flux are subtracted from the V and B 
band fluxes using the R band magnitude (ttir = 17.21) 



from |Urry et al. (2000) 
liptical galaxies (z 



and t he spectral template o f el- 
2) from |Fukugita et al. (1995)| The 
host-galaxy-subtracted fluxes are corrected for Galactic 
extinction using As = 0.058 mag and t he interstel- 
lar reddening curve by |Cardelli et al. (1989)[ updated by 
O'Donnell (1994)| The exposure ('J') and observed flux (F) 
in each optical band are shown in Table H] for each of the 
two nights. The contribution of the host galaxy to the other 
filters is negligible. 



4.2.4. NRT Results 



The NRT (jTheureau et al. 2007) . a meridian transit tele- 
scope with a main spherical mirror of 300 m x 35 m, 
measured the 11-cm flux from H 2356— 309 on 29 different 
days between June through December 2005 as part of an 
on-going monitoring program. The nearest observations to 
the XMM-Newton pointings were performed on the three 
nights of June 11 — 13, 2005. The observed 11-cm flux was 
quite low and not strongly variable during 2005, with an 
average of 7.5 ± 2.0 mJy. 



5. Discussion 

The broad-band SED of H 2356-309 is shown in Figure E 
and Figure [9] for the RXTE campaign in 2004 and for each 
of the two XMM-Newton observations in 2005. All of these 
SEDs have the same shape as most VHE HBL, indicat- 
ing a double-peaked structure with a roughly similar en- 
ergy output in both the keV and VHE range. It should 
be noted that the optical, X-ray and VHE measurements 
from each of the three campaigns (RXTE, XMM-Newton 
on MJD 53534 and XMM-Newton on MJD 53536) are very 
similar. The measured X-ray flux from 2 to 10 keV differs 
by only ±15% from the mean value, and the X-ray spec- 
tral slopes are consistent above 2 keV. The annual VHE 
fluxes differ by less than ±30% from the mean and the an- 
nual photon indexes are constant in time. As the VHE and 
X-ray spectral slopes are all consistent, and only modest 
changes in the VHE and X-ray fluxes are observed, it is ex- 
pected that modeling of each of the three SEDs will result 
in similar physical parameters. 

The aforementioned characterization of these three 
SEDs i s performed using a sim ple, time-independent SSC 
model ( jKatarzyfiski et al. 200Tj ). In this SSC scenario, the 
optical through VHE emission is generated by a spherical 
emission region of radius R relativistically propagating with 
a Doppler factor S with respect to the observer through 
a homogeneous magnetic field B. The emission region is 
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Fig. 8. SED of H 2356-309 in different epochs. The black triangles represent the SED from the 2004 RXTE campaign. 
The closed-blue and open-red circles represent the SED from the 2005 XMM-Newton observations on MJD 53534 and 
53536, res pectively. The light-gray tr iangles correspond to the BeppoSAX X-ray data during the high state observed in 
June 1998 (Costamante et al. 2001) In the VHE band, the black triangles and red circles are the time-average spectra 
observed in 2004 and 2005, respectively, and are not corrected for the effects of EBL absorption. The curves are fits to the 
respective SEDs of a single-zone homogeneous SSC model described in the text. The solid and dashed curves represent 
the same model with and without, respectively, the EBL effects included (Franccschini 2008). The SSC model fit to the 
2004 RXTE campaign data is identical to the fi t to the MJD 53536 data (red curve). H 2356-309 does not appear in the 
Fermi-LAT bright- AGN list|(Abdo et al. 2009)1 an d only poorly-constrained, non-simultaneous spectral data exist in the 
first Fermi-LAT catalog (Abdo et al. 2009)| The NRT (11 cm) measurement falls well below the minimum flux density 
of the figure and is not shown for the sake of clarity. For the luminosity axis, the measured redshift of z=0.165 is used 
with H = 70 kms" 1 Mpc" 1 , Cl M = 0.3. 



filled with relativistic electrons that generate the lower- 
energy peak at X-ray frequencies via synchrotron radiation 
and the higher-energy peak at VHE frequencies through 
inverse-Compton scattering of the synchrotron photons off 
the same electrons. The high-energy electron distribution 
between Lorentz factors j m i n and 7 max is assumed to have 
a broken power-law shape with normalization K and break 
7b. The energy index of the electrons is n\ below 7b, and 
ri2 above 7b. A low-energy cutoff in the electron distribu- 
tion of 7 m i n = 1000, is chosen to prevent an increase of the 
inverse-Compton emission in the MeV-GeV domain, that is 
not compatible with the non-contemporaneous upper limit 
on the blazar's MeV-GeV flux derived from the first 5.5 



months of Fermi-LAT data (Abdo et al. 2009) The major- 
ity of the radio emission is assumed to be produced from 
regions further out in the jet from the core. Therefore the 



NRT radio measurements are not used in the following SED 
modeling of the blazar, but can be considered as upper lim- 
its for the blazar's radio flux. 

As H 2356-309 has a redshift of z=0.165, its ob- 
served VHE spectrum is strongly affected by the absorp- 
tion of VHE photons on the EBL (see, e.g., the review 
of lAharonian 2001]) . These interactions (7 VHE 7 EBL — > e + e~; 
IGould fc Schreder 1967|) along the line-of-sight create an 
energy-dependent opacity t{z,E) that is imprinted on the 
observed VHE spectrum (F ohs (E) = F int (E) e~< z ' E )), To 
remove this absorption in the SED modeling the EBL den- 
sity of [Franccschini (2008) is used to calculate the opti- 
cal depths relevant for these observations of H 2356— 309. 
Fitting the de-absorbed time-averaged VHE spectrum to 
a power-law function (dN/dE ~ _E _rint ) yields a relatively 
hard intrinsic photon index of L; nt = 1.97 ± 0.15. Similar 
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Table 5. Results of SSC fits to the SED of H 2356-309. 



Log E [eV] 
11.6 11.1 





XMMi^ 


XMM 2 ° 


RXTE 1 " 


AH06^ 




lo 


1 8 

lo 


1 s 
lo 


lo 


B [G] 


0.16 


0.16 


0.16 


0.16 


R [cm] 


7.5 x 10 15 


7.5 x 10 15 


7.0 x 10 15 


3.4 x 10 15 


K [cm" 3 ] 


8.0 x 10 4 


8.0 x 10 4 


2.0 x 10 4 


1.2 x 10 4 


Tmin 


1.0 x 10 3 


1.0 x 10 3 


1.0 x 10 3 


1.0 x 10 3 


7b 


1.0 x 10 5 


1.0 x 10 5 


2.5 x 10 5 


2.5 x 10 s 


7max 


3.0 x 10 6 


3.0 x 10 6 


3.0 x 10 6 


3.0 x 10 6 


m 


2.3 


2.3 


2.0 


2.0 


112 


3.6 


3.5 


3.4 


4.0 



a The XMMi and XMM 2 results refer to the fits to the data 
from MJD 53534 and 53536, respectively. Using the XMM- 
Newton optical spectrum from MJD 53536 instead of the average 
ROTSE-IIIc flux from 2004 for the RXTE campaign yields the 
same SSC fit parameters as shown in XMM2. 

b This is the result of the original SSC modeling of 
the RXTE campaign SED of H 2356-309 published in 
Aharonian et al. 2006b~l The RXTE 1 " results use the same multi- 
wavelength data as AH06 along with an improved calibration of 
the VHE data. 



values of r;„t are found when fitting the de-absorbed data 
from the 2004 and 2005 VHE spectra used in the SSC mod- 
cling. It should be noted that the approximate slope of each 
of the de-absorbed SSC curves (dashed lines in Figures \E\ 
and [5]) is slightly steeper than the power-law fit to only the 
HESS EBL-de-absorbed data. This is due to the inclusion 
of all multi-wavelength data in the curved SSC fits, whereas 
the power-law function is fit to only the HESS data. 

The parameters of the SSC fit to the SED from each 
of the XMM-Newton campaigns are shown in Table [5J For 
the modeling of the RXTE campaign, the XMM-Newton 
optical spectrum from MJD 53536 is used instead of the 
average ROTSE-IIIc flux from the 2004 season, since they 
are consistent and the former more strongly constrains the 
fit. In addition, the ROSTE-IIIc R-band measurement is 
modeled from white light, whereas the XMM OM data are 
from filtered light and hence are more appropriate for SED 
modeling. The data from the RXTE campaign are well de- 
scribed by the parameters of the SSC fit to the data from 
the second XMM-Newton night (MJD 53536) and modi- 
fication of these parameters does not yield an improved 
X 2 . Indeed the fit parameters do not vary greatly between 
the three campaigns, as expected from the similarity of the 
measured quantities. This may suggest that a steady state 
of the blazar was observed, and that this state can be ex- 
plained by an SSC scenario. The slight hardening of the 
XMM-Newton spectrum between MJD 53534 and 53536 is 
easily accounted for by only a minor variation in the high- 
energy index of the electron distribution. 

It should also be noted that the parameters of the 
SSC fits indicate a cooling time scale of ~1.5 x 10 5 s 
(~4 h), which is comparable to the source crossing time 
(t esc = R/c = 2.5 x 10 5 s), but is considerably shorter than 
the monthly time scale of the weak variations observed in 
the VHE band. However, VHE flux variations on time scales 
comparable to the cooling time (e.g. daily), and of approx- 
imately the same magnitude, cannot be ruled out due to 
the low average VHE flux. 




26 26. 2 

Log i/[Hz] 



Fig. 9. SED of H 2356-309 zoomed in on the VHE band. 
The data and model curves are identical to those in 
Figure [5] 



The comparison of the results of the fits to the RXTE 
campaign m 2004 (see XMM 2 in Tab le [5]) to the previously 
published results (jAharonian et al. 2006b[) shows some dif- 
ferences. These differences are not primarily due to the 
renormalization (^50% increase) of HESS flux from this 
epoch. Rather they are largely due to the improved knowl- 
edge of the optical spectrum of the blazar from the XMM- 
Newton OM results. Using the ROTSE-IIIc flux instead of 
the XMM-Newton OM spectrum results in a different set 
of fit parameters (see RXTET in Table [5]) , very similar to 
the previous modeling (see AH06 in Table[S|). However, this 
fit (RXTE t ) falls significantly below the measured XMM- 
Newton optical spectrum. As the optical portion of the SED 
is not expected to change dramatically, the characterization 
presented here (XMM 2 ) should be used as a reference for 
future studies. 



6. Conclusion 

Long term HESS observations of H 2356— 309 have clearly 
confirmed the initial VHE detection from 2004. The source 
has a relatively low VHE flux (1.6% Crab) that is weakly 
variable (factor of ^2) on time scales of at least months. 
Although the flux is variable, no changes in the VHE spec- 
tral slope (r w 3.1) are found. After accounting for the 
effects of the absorption of VHE photons on the extra- 
galactic background light, the intrinsic VHE spectrum of 
H 2356— 309 is found to be hard (Fi nt ~ 2.0). Interestingly, 
the X-ray spectrum is also very hard, with the synchrotron 
peak located at ~1 keV, despite the source being in a his- 
torically low state in X-rays. While the synchrotron peak of 
H 2356—309 is located at higher than typical energies and 
the X-ray and VHE spectra are comparatively hard, the 
SED can be reasonably modeled using a simple one-zone 
SSC scenario. This could be expected since the luminos- 
ity of both SED peaks are similar, and their separation in 
frequency is not extreme. 
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